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Abstract. Reaction— diffusion equations with a fractional Laplacian are reduced 
near a long wave Hopf bifurcation. The obtained amplitude equation is shown to 
be the complex Ginzburg-Landau equation with a fractional Laplacian. Some of 
the properties of the normal complex Ginzburg-Landau equation are generalised 
for the fractional analogue. In particular, an analogue of Kuramoto-Sivashinsky 
equation is derived. 
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1. Introduction 

Random processes characterised by Levy flights have been discovered several decades 
ago. Since then similar processes have been observed in numerous natural phenomena: 
diffusion and advection in fluids [T], in particular turbulent flows [2] and wave 
turbulence [3], motion of animals [4], balance control in humans [5] and even progress 
of seismic foci [6] . At the macroscopic limit Levy flights are modelled by a fractional 
Laplacian operator. General properties of such processes are reviewed in [7]. Often 
non-linear kinetics, such as ( but not only ) chemical reactions create an intricate 
interaction with the diffusion process, especially at an instability threshold ^-|10j. 
Thus an equation combining a fractional diffusion operator ( a Laplacian fraction 7/2 
with 1 < 7 < 2 ) and non-linear kinetics is the simplest model to capture the basic 
effects of such an interaction. 

Understanding of pattern emergence and formation in normal reaction-diffusion 
systems near a Hopf bifurcation point was achieved by means of amplitude and phase 
diffusion equations |ll|-|14j. The current work follows the course of reduction of the 
fractional reaction-diffusion model near such a bifurcation in order to obtain and study 
the fractional analogues of complex Ginzburg-Landau ( amplitude ) and Kuramoto- 
Sivashinsky ( phase diffusion ) equations. 
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2. Reduction near bifurcation point 



Consider a two species fractional reaction-diffusion system 

-f(n), (1) 



dn / di \ / ®N"i 



dt V d2 J \ J)7.^^ 

where the Laplacian fractional counterpart is of a generally distinct order for each 
species: 

and n,f and dj are the species concentration vector, kinetics function and diffusion 
coefficients, correpondingly. Suppose there exists a uniform steady state ng satisfying 
f(no) = 0, in whose close vicinity f varies according to a sensitivity matrix (Vf)jfc = 
dfj/duk, j,k S {1,2}. Then Hopf bifurcation occurs when its trace vanishes. Let us 
split the matrix as 

Vf = Vfo + Q ^ , trVfo - 0, e< 1, < 0(1), (3) 

rescale 

n{x,t)^-N{^,to,t2,...;e), ^ ^ Sx, tj ^ eH, j ^ 0,2, . . . (4) 

and expand 

oo 

N r-.no + ^6jNj{£„to,t2,...), Sj^Sj{e). (5) 

For the leading order reduction only the first slow temporal and spatial scales will 
be used. Substitution ^ into ([ij and scrutiny of the resulting system lead to the 
following choice of 6 and Sj. If "fi — j2 — 7, the scales are 6j — and S'' — e^, where 
the latter ensues by 

Dl^y{x)=6^^l^y{^/S). (6) 

Then at order 0{Si) a system of linear homogeneous equations at the bifurcation point 
is obtained: 

dNi 



„^ VfoNi = 0. (7) 

dto 



Ni = A($,i2)e'^*''vi+c.c., X = iuj, (8) 
1 \ ^ itj-V/ii 2 



Its solution is 



v,= ^^ j, "=^^, .^^detVfo. 

At subsequent orders the system is not homogeneous. At order 0(63) secular non- 
homogeneous terms coerce a solvability condition for A, alias fracitonal amplitude or 
Ginzburg-Landau equation: 

OA _ t^_^ , fd2 + di , .d2 - di Vf: 



dt2 2' 



['^^^'^^-^)%,A-sA\A\^ (9) 
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where s is constant and depends on lo and derivatives of f up to third order at the 
bifurcation point. To obtain the normal form rescale the variables as 

d2 - di V/u ^ 3s 



d2 + di Lu 5Rs 
8 A 

— ^A+{1 + m)S)f^| A - sign(3fis)(l + f3i)A\A\\ (10) 

If however 71 < 72, the expansion scales 6j remain the same up to third order, where 
the anomalous term first appears: 6j — , j € {1,2,3}. The spatial scale choice is 
according to the activator exponent 6'''^ — . Then 1 + 272 /71 > 3 and the anomalous 
term of the inhibitor will be neglected at order 63. Actually, it will appear only at 
order k + 1, where k is the greatest integer satisfying 1 + 272/71 > k. The amplitude 



equation coincides with (10 1 upon setting 7 = 71 and c?2 = 0. When 71 > 72, the 
same holds with 7 = 72 and di = 0. Thus the fractional analogue of the complex 
Ginzburg-Landau equation is 
dA 

— = A+(l+m)^;;|v4-sign(3f^s)(l+/3^)A|A|^ 7 = min{7i, 72}.(11) 

This equation was formerly derived in 15J in the problem of nonlinear oscillators' 
dynamics with long range interactions. In the current paper the properties of the 
super-critical version ( 5Rs > ) are studied. 

3. Symmetry properties 

The normal complex Ginzburg-Landau equation possesses translational space x — > 
X + Cx, time r ^ r + and phase A ^ A exp(ic) symmetry. Moreover, solutions 
within the class of modulated waves 

^(a;,T) = B(r)e^(9^-"='^', r ^ x - vt, q,m,veR (12) 

are connected by a similarity transformation within the family (a— (3) /{l+af3) = const 



[TT] . With the introduction of anomaly (11) loses the Galilean invariance, and the 



symmetry within the class of modulated waves is preserved only for w = 0. A solution 
with given (a,/3) connects to a solution with {a' , /3') hy B = aB' , r — br' . Then 
q' = bq, 

1 + a' 13' 1 + 



1-1- a/3 l-fa' 
1 -I- a'^ 

1 -t- a2 ^ ' 

b-< = ^ + + (13) 
4. Variational formulation 

For the special case a — (3 the normal complex Ginzburg-Landau equation can be 
obtained by variation of a functional [TT]. Rotation A 1-^ Aexp(— z/3r) gives 
f)A 

— ^{l + ^f3){A + Dl^A-\A\'A). (14) 
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Then the functional must be adjusted to yield the fractional Laplacian operator: 

U{x, T)dx, 

ll|2 I 1|^|4 Bec(.7/2) [dA* d r A{OdQ 



1 



dA* sign(a: — C) 



rfC + c.c. \ + c 



(15) 



It is possible to choose the constant c so that T converges, as the function A is 7-fold 
differentiable. The first variation is 



Since 



(5T = - 
dA _ 



/OO 
[A- A^A* + 'B'I^^A)5A*dx- 
-00 



- c.c. 



SU 



all solutions of ( 14 1 decay asymptotically in time 

dU 



dT 



dT 



-dx 



2 


OA 


1 + J -00 


dt 



dx < 0. 



(16) 
(17) 

(18) 



5. Stability of traveling waves 



Traveling waves Aq{x,T) — pexp{i{qx — wt)) are a sub-class of (12) and comprise an 
important solution family of (11 1. By 



^7 piqx 



it is straightforward to show that 

p' = l-\q\\ ^ = (3-{f3-a)\q\''. 

5.1. Spatially homogeneous oscillations 



Linearising (111 about Aq, i.e. 

^ = e^'^"'(l + U + w), u,weM, |m|,|u|<1, 
splitting into a real system and neglecting non-linear terms of u and u, 



d_ 



1 

(3 



1 

a 



-a 
1 



The eigenvalues A of a normal disturbance 

Xr+ikx 



Ul 



satisfy 

+ 2A(1 + |A:|T) + |fcp ((1 + a2)|fc|7 + 2(1 + a/3)) = 0. 
Since Ai + A2 = —2(1 + |fc|''') < 0, the disturbance is unstable if 
A1A2 = \k\-' ((1 + a2)|fc|T + 2(1 + a/3)) < 0, 

which gives a set of unstable wave numbers when 1 -t- a(3 < 0: 

^ l + a/3 



< l/cl < fc„ 



(19) 
(20) 

(21) 
(22) 

(23) 

(24) 
(25) 

(26) 



1 + a2 ■ 

Thus the instability domain in the (a, /3) plane coincides with the normal Benjamin 
Feir domain. 
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5.2. Arbitrary wave 

The evolution of a small perturbation a(a;, r) about an arbitrary wave solution Aq is 
governed by 

1^ = a + (1 + za)D^^|a - (1 + tf3){2a\Aq\' + a*Al). (27) 

Without loss of generality the base wave may be taken one-dimensional. However the 
disturbance should combine a longitudinal and transverse waves: 

a = Ag+fc(T)e*(«+'=-)^+''=''^ + yl,_fc(T)e'(«-'=-)^-*'=''^. (28) 
The resulting system of equations is 

wherein 2tg is a 2 x 2 matrix whose entries depend on and ky. With Aq±k = 
A± exp((A =F jcc)r) the eigenvalues must satisfy a quadratic equation 
/ X~iw- -2ti2 \ 
det = 0. (30) 

\ -21^2 X + iw-^22 J 

Below some particular relations between the longitudinal and tranverse disturbance 
wave numbers are considered, and (30 1 is expanded appropriately at the corresponding 
limits. 

Note that the underlying wave Aq is neutrally stable: for k^ — ky ^ Q the 
eigenvalues are Ai = 0, A2 = —2(1— \q\^) < 0, thus rendering the long wave 
disturbances of special interest. For small ratios kx/q, ky/q one can distinguish 
between two qualitatively distinct cases: 

(z) 0(%/g) ^ 0(fc,/g) ^ 0(i.) 

(zz) 0{ki/q) ~ 0(fc2/g2) ^ o(^) ' « i^^i 

First suppose that ( [3T| ) holds. Then the expansion is taken up to order 0{v^) because 
at 0{v) the real part of Ai vanishes: 

5RAi ^l\qV (-{I + a(3) L - 1)| + M^j + (1 + P')^^'^ ) + (32) 



Thus if 1 + a(3 < 0, instability is immediate for any wave (20 1. For 1 + a(3 > solving 
to leading order the inequahty 3fiAi > yields 

wherein coercing positiveness of the right-hand side gives a subset of unstable wave 
numbers 

q,n<\q\<l, q-: = l+ ^.l^^l - (34) 

7 — 1 1 -f ap 

Then in the (A:^, fc,,) plane the instability region is bounded by two intersecting straight 
lines. For traveling waves within this subset pure longitudinal disturbances ( fc^ = ) 
have the highest growth rate, recovering the Eckhaus instability criterion for the 
normal Ginzburg-Landau equation. 

Now suppose (31 ;i) holds. In this case an expansion to order O(i^) suffices: 

^Xi^-^\qp{l + a[3)^ + 0{i^^). (35) 
2 (j^ 

Here the instability ensues only for 1 -|- aP < 0. 
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6. Phase diffusion equation 

At the opposite limit of small ratios q/k^, l/ky no new instability criteria emerge. The 
spatially homogeneous oscillation Aq is unstable within the same region 1 + aj3 < 
with respect to disturbances ( 26 ) . The evolution of perturbations near the domain 
boundary is described by a fractional non- linear phase diffusion equation ( Kuramoto- 
Sivashinsky equation fractional analogue ). Define < e <C 1 so that 1 + af3 — —e. 
By (26) the spatial coordinate scale is x = e^/'^x. To find the appropriate temporal 
scale take |fc|^ = Ke and expand (24) in powers of e. The resulting approximation is 

1 



Ai 



K 



Hence the temporal scale is T2 = e^r. Using ^ and rewriting (111 with x and T2 



Taking 



where 



. A + eil--il + e))^l^A-{l + zmA\'A. 



'•{X,'r2) = 1 + rj{x,T2), 93 = ^ eVj(x, ^2), 



substituting into (37 1, dividing by exp{iip) and using the expansions 

2 

00 / 00 \ 



00 ^ / QO \ 



(36) 

(37) 
(38) 
(39) 

(40) 



it is possible to collect powers of e. At order O(e^) the following equation is obtained: 



dipi 
9t2 



i)~2)i^,l^i.(41) 



The operator {D7 i)^ is defined in Fourier space by 



\x\' 
(D 



(42) 



and cannot be related in a simple way to the operator S^Jp because the order 



27 exceeds the range of definition of S''' 



lx|- 



Note that the coefficients of the linear 



terms D^^^cpi and (S)|'^|)^(/ji are consistent with (36 1). To bring (41) into parameter 
independent form define 

T = T2T, x^bx, ip = aipi. 

Then 



(43) 



dip 
b'< ■ 



2/32 



a = P + 



1 



(44) 



Except for replacing the Laplacian by fractional operators the resemblance to the 
Kuramoto-Sivashinsky non-linear phase diffusion equation is obvious. 
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7. Numerical simulations 



The fractional complex Ginzburg-Landau equation (11) and fractional Kuramoto- 



Sivashinsky equation (44 1 have been solved numerically by a pseudo-spectral method 
with time integration in Fourier space, Crank- Nicolson scheme for linear operators 
and Adams-Bashforth scheme for nonlinear ones. Periodic boundary conditions, and 
if not specified otherwise, small amplitude random data as an initial condition have 
been used. Due to the extremely diverse dynamics of both CGL and KS equations 
[TTl [T7] the current section is limited to the most interesting regimes that show the 
difference between the normal and fractional equations. 

In one spatial dimension phase and amplitude turbulence, two regimes of complex 
behaviour known for the normal COLE [llj . are comparable with the anologous 
dynamics of (11). Figure [T| shows phase turbulence regimes for different values of 
7 with a and (3 close to the Benjamin- Feir instability threshold. As initial condition 
a spatially homogeneous state ReA — const, ImA = const, \A\ = 1 was used with 
small amplitude random noise added. The upper subfigure corresponds to 7 = 2.0, i.e. 
normal CGL equation with its typical picture of phase turbulence and small amplitude 
modulations near \A\ — 1, exhibiting spatio-temporal chaos in the form of splitting and 
merging "cells" [11]. This type of dynamics is described by the Kuramoto-Sivashinsky 
equation ( see below ). The middle subfigure shows a phase turbulence regime for 
7 = 1.7. In this case, together with merging and splitting cells, long-living traveling 
shocks form, propagating in different directions. These shocks "absorb" the adjacent 
"cells". The lower subfigure corresponds to 7 = 1.5. Here a single shock forms in the 
whole domain and travels with a constant speed. The shock tails exhibit weak chaotic 
modulations. As shown below, all above regimes can be captured by a fractional 
KS equation, describing phase dynamics near the Bcnjamin-Feir instability threshold. 
With increasing distance from the threshold the shocks self-accelerate and trigger a 
transition to an amplitude turbulence regime, as shown in figure |2] 

Figure [3] shows amplitude turbulence regimes for different values of 7. The 
usual amplitude turbulence dynamics comprising traveling hole solutions with a weak 
component of the phase turbulence is observed for 7 = 2 ( upper subfigure; see also 
|llj ). As 7 decreases, the phase turbulence component grows stronger, and for 7 1.0 
( lower subfigure ) a combined phase-amplitude turbulent regime is formed with no 
traveling hole solutions. 

Experimental observation of phase and amplitude turbulence regimes described by 
a ID FCGL equation require special arrangements in order to make the experimental 
system effectively one-dimensional. A more generic and easier to study in experiments 
are two-dimensional systems, described by a 2D FCGL equation. Probably the most 
remarkable solution of a normal CGL equation in 2D is a spiral wave ;11]. Therefore 
the effect of anomalous diffusion on spiral wave dynamics is of prime interest. 

Numerical simulations of 2D FCGL equation have been performed for parameter 
values, corresponding to the formation of spiral waves in normal CGLE, by means of a 
similar 2D pseudo-spectral code. Figure |4] shows snapshots of the solutions for different 
values of 7. Figures [4]2, c? correspond to 7 = 1.9. A single spiral wave is formed in 
the whole domain, akin to the normal case. In figures |4|), e ( 7 = 1.8 ) spiral-like 
defects form, where the core of each defect occupies a certain domain with domain 
walls between them. However, the domain walls in this case are partially dissolved 
and the defects ( spiral cores ) move in a chaotic manner. With further decrease of 
7 ( figures [4]:, /, 7 ~ 1.05 ) the number of defects decreases and the domain walls 
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between them are almost completely dissolved. The system exhibits several domains 
with almost spatially homogeneous oscillations with different phases. 

Finally, the fractional Kuramoto-Sivashinsky equation (44 1, describing phase 
dynamics near the Benjamin-Feir instability threhsold, has been solved by a pseudo- 
spectral code with periodic boundary conditions and small amplitude random data as 
an initial condition. Figure [5] shows spatio-temporal solution diagrams for different 
values of 7. Figure [sji, corresponding to normal diffusion {j — 2.0 ), shows chaotic 
spatio-temporal dynamics of merging and splitting "cells", as in phase turbulence of 
the normal CGLE shown in figure [ij upper ). Figures [sj), c conform to 7 = 1.7 and 
7 = 1.6, respectvely. Along with merging and splitting cells traveling shocks appear, 
absorbing and emitting cells. As 7 decreases, the shocks become more frequent and 
pronounced, and propagate faster, similarly to figure [l]( centre ). When 7 decreases 
below a certain, domain length dependent threshold, a single traveling shock is formed 
in the whole domain, as exemplified in figure |6] a spatio-temporal diagram (a) and 
a few solutions in successive moments of time (6). The shock travels with a constant 
speed while its tails exhibit weak, chaotic spatio-temporal modulations. A similar 
shock formation is seen in figure [l] ( lower ) . 

Note that the shock amplitude grows both with the decrease of 7 and increase 
of the computational domain. Below some threshold of 7 the shock accelerates with 
its amplitude growing exponentially. An asymptotic analysis revealed large amplitude 
asymptotic solutions of ( 44 1 in the form 

^ = a(T)/(x-C(r)), (45) 

where / is an odd periodic function, a(r) grows exponentially, and the instant speed 
d(/dT is proportional to a(T), as confirmed by numerical simulations. 




Figure 1. Amplitude \A\ in a phase turbulence regime as described by FCGL: 
spatio-temporal diagrams of numerical solutions of for 7 = 2.0, /3 = 1.2 
( upper ), 7 = 1.7, /3 = 1.1 ( centre ), 7 = 1.5, fS = 1.05 ( lower ); a = -1.0. 
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Figure 3. Amplitude A in amplitude turbulence regi me as described by FCGLE; 
spatio-temporal diagrams of numerical solutions of for 7 = 2.0 ( upper ), 
7 = 1.5 ( centre ), and 7 = 1.0 ( lower ); a = 1.0, /3 = —1.3. 



8. Discussion and conclusions 



A reaction-difTusion system governed by Levy flights has been reduced near a long 
wave ( Hopf ) bifurcation point. The reduced system is described by a complex 
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a) b) c) 




' o' o' ' 

250 Q 250 250 



Figure 4. Solution snapsliots of 2D FCGLE: 5RA ( upper ) and |^| ( lower ) for 
a = 1.5, (3 = -0.6, and 7 = 1.9 (a),(d); 7 = 1.8 (b),(e); 7 = 1.05 (c),(f). 

Ginzburg-Landau (FCGL) equation with a fractional order Laplacian. 

The fractional analogue does not inherit the Galilean invariance, a property well 
known for a normal Ginzburg-Landau equation, yet the similarity of modulated wave 
solutions along the family of curves (a — /3)/(l + a(3) = const is retained. Another 
property common to normal and anomalous diffusion is the variational formulation in 
the special case a — (3. 




^ ID00 



Figure 5. Spatio-temporal dynamics of FKS equation l |44[ l for 7 = 2.0 (a), 
7 = 1.7 (b), 7 = 1.6 (c). 
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Figure 6. Solutions of FKS equation l |44| for 7 = 1.5: spatio-temporal diagram 
(a) and solutions at sucessive moments of time (6). 



Similarly to a normal complex Ginzburg-Landau equation, the fractional analogue 
possesses the family of solutions in the form of plane traveling waves. For a uniform 
oscillation the instability domain in the a-P space coincides with the normal Benjamin- 
Feir domain. The instability with respect to two-dimensional disturbances ensues for 
the whole unit circle within this domain and for a 7-dependent subset otherwise. 

Near the Benjamin-Feir domain boundary the system dynamics is described by a 



fractional analogue of the Kuramoto-Sivashinsky equation ( 44 1 for the phase evolution 



Numerical solution of one-dimensional FCGL equation for decreasing values of 
Levy exponent 7 reveal the appearance of traveling shock waves in the regime of 
phase turbulence, whereas amplitude turbulence exhibit a stronger phase turbulence 
component. The decrease of 7 in two-dimensional solutions leads to destruction of 
spiral waves and formation of defect chaos. 

Solutions of FKS equation for decreasing 7 exhibit a transformation of chaotic 
dynamics of merging and splitting cells, typical of systems with normal diffusion, to 
traveling shocks. The same transition is observed for the FCGL equation. When 
7 decreases below a certain threshold, the shocks self-accelerate, their amplitude 
grows exponentially and the solution blows up. This phenomenon corresponds to the 
transition to amplitude turbulence and has been confirmed by numerical simulations. 

Some remarks on the possible ways to control the anomalous diffusion effects in 
experiments are in order. For a reaction-diffusion system set up in a liquid layer with 
turbulent mixing the means to control the Levy exponent, characterising turbulent 
diffusion, would be to control the mixing intensity and possibly the type of turbulent 
flow. For a reaction-diffusion system on a catalytic surface, in which the surface 
super-diffusion would be a consequence of the possibility for molecules to make long 
jumps over the surface through the gas phase, the Levy exponent might be controlled 
by temperature change or control of the gas flow neax the surface by affecting the 
adsorption bonds, say, by light irradiation. 
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